Entanglement and symmetry effects in the transition to the Schrodinger cat regime 
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We study two-spin entanglement and order parameter fluctuations as a function of the system size 
in the XY model in a transverse field and in the isotropic XXX model. Both models are characterized 
by the occurrence of ground state degeneracy also when systems of finite size are considered. This 
is always true for the XXX model, but only at the factorizing field for the XY model. We study 
the size dependence of symmetric states, which, in the presence of degeneracy, can be expanded 
as a linear combination of broken symmetry states. We show that, while the XY model looses its 
quantum superposition content exponentially with the size N, a decrease of the order of 1/N is 
observed when the XXX model is considered. The emergence of two qualitatively different regimes 
is directly related to the difference in the symmetry of the models. 
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I. INTRODUCTION 



The transition between the microscopic and macroscopic worlds is a fundamental issue in quantum theory both 
from the point of view of foundations of physics and of the application to quantum computation [l|, |2[ . 

Spontaneous symmetry breaking (SSB) indicates a situation where, given a symmetry of the Hamiltonian, there are 
eigenstates which are not invariant under the action of this symmetry, unless a term is added which explicitly breaks 
the symmetry. Usually, when the control parameter reaches a critical value, the lowest energy eigenstate keeping the 
Hamiltonian symmetry is no longer stable in the presence of infinitely small perturbations, and new stable solutions 
appear which are not symmetric. SSB leads naturally to a degenerate manifold of ground states. 

Symmetry breaking usually occurs in the thermodynamic limit, when superselection destroys quantum coherence. 
Important exceptions are the XXX Heisenberg model and the XY model in a transverse magnetic field at the particular 
value of the field where ground state factorization occurs @. In these cases, ground state degeneracy occurs for any 
size of the system, and it is therefore possible to explicitly study the transition from quantum to classical behavior. 
How entanglement is affected, in the thermodynamic limit, by the presence of a term which explicitly breaks the 
symmetry, ha been discussed by Syljuasen Q and by Osterloh et al. 5]. Here, we face the problem from another 
point of view, starting from small systems, and then increasing the size until the thermodynamic limit is reached. 

The XY model in a transverse field has been introduced in the early sixties and solved by Katsura @, by means 
of the Jordan- Wigner transformation, formerly introduced by Lieb, Schultz, and Mattis [7[. Subsequently, the corre- 
lation functions were investigated in great detail by Barouch and MacCoy 8], who found the existence of a second 
critical value of the transverse field separating qualitatively different behaviors of the correlation functions. Later 
on, Kurmann, Thomas and Muller |3j discovered ground state factorization for a large class of spin models. In the 
particular case of the XY model, the field at which factorization occurs is exactly the critical field of Barouch and 
MacCoy. Recent interest has been devoted to the study of entanglement properties of many-body systems undergoing 
a quantum phase transition [jj-[i~2|. As shown in Ref. 13 1, the critical point turns out to separate two regions with 
qualitatively different bipartite entanglement. It has been shown in Ref. 14] that, in the vicinity of the factoriz- 
ing field, the range of concurrence diverges, and that such divergence corresponds the appearance of a characteristic 
length scale in the system. Recently, the conditions for the existence of the factorizing field for models with long-range 
interaction have been determined 15 1, and the study of this effect has been extended to dimerized chains (lfjj . 
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The two systems we wish to investigate belong to different classes of symmetry. While the XXX Heisenberg model 
has the SU(2) continuous symmetry, i.e. the Hamiltonian commutes with the total spin along any possible direction, 
the XY model is invariant under parity transformations and possesses the discrete symmetry. 

It is commonly accepted that purely quantum effects arc not observable on the macroscopic scale, except for super- 
conductivity, superfluidity. On the other hand, quantifying entanglement (perhaps the most genuine manifestation of 
quantum properties) as a function of the system size represents a fundamental issue [13]. Here, we wish to investigate 
in detail the two-spin entanglement dependence on the total number of spins for these models. In particular, we shall 
derive the difference in the size effects due to the difference in the system symmetry. To be more specific, in the case of 
a discrete symmetry there is an exponential entanglement decrease, while with a continuous symmetry entanglement 
shrinks linearly with the growth of the system size. 

Since the XY model has been studied through the last four decades, and results are scattered over a vast literature, 
for convenience we shall give here a brief survey of the main results, with the primary aim of focusing on the existence 
of the factorizing field and its independence from the system size. 

The paper is organized as follows. In Sec. |n] the XY model in a transverse field is discussed. We give special 
emphasis to the finite size solution with the scope of enlightening the emergence of the factorizing field as a size- 
independent degeneracy point. Furthermore, by means of the finite size picture, we are able to explain in a simple 
way the appearance of spontaneous symmetry breaking in the thermodynamic limit. In Sec. IIIII we describe briefly 
the structure of the ground state for the isotropic Heisenberg (XXX) model. Even if, for any finite number of spins 
the ground state manifold has finite dimension, an over-complete set of states can be introduced that allows to study 
the microscopic-to-macroscopic transition. In Sec. II VI we derive the value of the concurrence for pairs of spins and the 
order parameter fluctuation in a superposition state as a function of the size systems both for the XY model and the 
XXX model. Finally, in Sec. |Vj results are discussed. In particular, we will focus on the influence of the symmetry 
in the different behaviors. 



II. XY MODEL 



Let us consider a chain of N spins 



H = /Z 



J 



(1+7) 



'i u i+i 



J 



(1-7) 



V V 

a l a l+l 



h<7? 



(1) 



where ct c arc the three Pauli matrices (e — x, y, z), and periodic boundary conditions (c^r +1 = a\) are assumed. In the 
following we will assume J = — 1 (ferromagnetic coupling). The above Hamiltonian is invariant under the Zi group 
of the rotations by tt about the z axis, since it commutes with the parity operator P = \\ l <jf. Due to this symmetry, 
eigenstates are classified depending on parity eigenvalue. This system is known to undergo a quantum phase transition 
at the critical point h c = 1. Below this value, in the thermodynamic limit, spontaneous magnetization along the x 
axis appears. 

Since the work of Ref. [7], the Jordan- Wigner transformation, defined through af = 1 — 2cjcj, er ; + = 
Y[j<i (l — 2c;C/J a, <7 ; ~ = Yij<i fl — ^4 c ij 4 i s introduced to map spins in spinless fermions. The transformed 
Hamiltonian is H = Hq — PH\ with 



N—l 



H = [(4 c i+i - c i4+i) + 7 (44+1 - c * c *+i) ~h(l- 2cJqJ 



Hi = 



[(4 



JyCi - c N c[ +7 



(cUt 



C N Ci 



(2) 
(3) 



Since [H,P] = 0, all eigenstates of H have definite parity, and we can proceed to a separate diagonalization of H 
in the two subspaces labelled by to P = ±1. Then, the complete set of eigenvectors of H will be given by the even 
eigenstates of H + = Hq — Hi and the odd eigenstates of H~ = Ho + Hi. Both for H + and H~ the diagonalization 
is obtained by first carrying out the space Fourier transform 



Ck 



—y 



e l « M Q, 



(4) 



where k — 0,1, ... , N— 1 in H , and k = 1/2, 3/2, . . . , N— 1/2 in H + , and then making the Bogoliubov transformation 



c fe = cos7? fc 77fe + zsin$ fe 7/I fe , 



(5) 
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with §k — —fi-k- Here, rj^k stands for rjN-k- The diagonalization condition implies for iDk 



7 sin A; ,„, 

tan 20* = —=-* -. 6 

h — cos k 



Eventually, we end up with the quasi-particle Hamiltonians 

N-l/2 



H 



+ = Afc Mm. - \ ) , (7) 

fc=l/2 ^ ' 



where the eigenvalues are given by 



A fc = 2^ ^/i-cos— fcj + 7 2 sin 2 — k. (9) 

The ground states of H + and _ff ~~ are the corresponding vacuum states with eigenvalues Eq = — J^k^if-z ^ k an< ^ 

The vacuum in the generic k mode is determined by r\k |0 ) = 0. While for every k ^ the Bogoliubov vacuum 
corresponds to an even state (the absence of quasi-particles implies zero or two particles), the mode k — plays a 
special role. In fact, the correspondent Bogoliubov transformation reads 

1 / h- 1 \ i ( h- 1 \ + 

^H^i^h+H 1 -^)^ (10) 

with the important consequence that the quasi-particle vacuum corresponds to a zero-particle state for h < 1 and to 
one-particle state for h > 1. The presence or the absence of the particle in the k = mode changes the parity of 
the state. Thus, for h > 1, the vacuum of H~ , because of its symmetry, does not belong to the set of eigenstates of 
H, while for h < 1 it becomes an eigenstate of physical interest. Above the h = 1, the odd state of lowest energy is 
obtained by adding one quasi-particle corresponding to the bottom of the energy band with energy A m j n = 2 (h — 1). 
This energy gap prevents the degeneracy even in the thermodynamic limit. 

A. Quantum phase transition and ground state factorization 

The change of symmetry of the vacuum of H~ is the very cause of the phase transition in the thermodynamic limit. 
Indeed, on the macroscopic scale the sum over k becomes an integral yielding Eq = Eq . Then, below the critical 
point he = 1 the odd and the even lowest eigenstates are degenerate, and the Hamiltonian symmetry is spontaneously 
broken, while, for h > he, due to the existence of the energy gap A m j n , the ground state keeps its parity (even). For 
h < he, because of superselection rules, the system is necessarily found in symmetry-broken states. 

As pointed out in Ref. [§}, below the critical point there are two different regions where two-body correlation 
functions can decrease monotonically or oscillate as a function of the spin distance, depending on the Hamiltonian 
parameters. These regions are separate, in the {h,~/} diagram, by the set of points satisfying h\ + 7 2 = 1. More 
recently, it has been shown that on this border line the ground state factorizes [23[, i.e. it can be written as j^^) = 

®i VpX with = (cost* |t*) ± sina||i)), where cos 2a = [(1 - 7) / (1 + j)} 1/2 . 

The existence of the factorizing field, originally derived by requiring only size-independent degeneracy between 
the lowest odd and even eigenvalues [3J, can be studied within the general solution of the model. By analyzing 
lowest odd and even eigenvalues of H in the symmetry broken region for finite N as a function of the transverse 
field, we observe a series of A^/2 level crossings for h — hi (see Fig. [I}. In correspondence of each hi the ground 
state changes its symmetry. The existence of such points has been discussed in Ref. [l8| and more recently in Ref. 
(l9| . and is responsible for the magnetization jumps reported in Ref. [2(J. In the thermodynamic limit, this kind of 
structure implies two different symmetry breaking mechanisms. For < h < hp, as N — > 00, the set {hi} of the 
degeneracy points becomes a denumerable infinity, while for hp < h < 1 there is a the usual symmetry breaking due 
to the vanishing of the gap. An interesting problem would be to check whether this is the microscopic mechanism 
responsible for the qualitative change in the behavior of the correlation functions above and below hp- 
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FIG. 1: Difference of energy between the lowest odd and even eigenvalues of H as a function of the transverse field for a 8-spin 
chain. Given the anisotropy amplitude 7 = 0.6, the factorizing field is /if = 0.8. As predicted, we observe N/2 level crossing 
point, the last of them being /if- 




While spontaneous symmetry breaking arises only for N — > 00, it can be seen from the previous analysis (see also 
Ref. [l9|]) that at the factorizing point hp degeneracy appears for any N. It is simple to show that Eq (hp) = Eq (hp) 
holds for any N, while the positions of all the other level crossing points hi change with N. Then, at the special field 
hp, the Hamiltonian symmetry is broken independently of the system size, and any linear superposition of the two 
symmetric eigenstates |a ± ) (\a + ) for the even eigenstate and \a~) for the odd eigenstate), is a possible eigenstate. 
Obviously, each one of the factorized states can be expressed as a linear combination of the two symmetric eigenstates 

=u+ \a + )±u- \a~), (11) 

1 /2 

with u± = [(l + cos Ar 2a) /2] . Notice that, for finite size systems, \^p) and \^fp) are not orthogonal, while 
{^'fI^&f} = in the thermodynamic limit. 



III. XXX MODEL 



The homogeneous (ferromagnetic) Heisenberg model is defined by the Hamiltonian 



JV-l 

Hxxx = - J E (°f °f+i + °f °f+i + *f+0 . (12) 
1=0 



with the boundary condition <j e N = <Tq. The model has been solved using the Bethe ansatz [21] . As far as the 
ground state properties are concerned, a simple argument can be introduced to show that, for any number of 
spins, any factorized state |$ (6, 0)) = ®/ [cos0|t) + exp (i<p) sm8 |T)] minimizes the energy. Given the invariance 
of Hxxx with respect to rotations of arbitrary amplitude (3 around any direction h 1Z (ft, h) = \\ l exp [ifiai ■ h], due 
to [Hxxx,Ti(P,n)] = 0, we shall restrict the attention on the particular state |<& (0, 0)) = |t,ti---)t)- It can be 
immediately seen that |$ (0, 0)) belongs to the ground state subspace of any of two-body terms of Hxxx, and then its 
energy represents the minimum achievable value. To make a link with the XY model, we could say that factorization 
point for the Heisenberg model corresponds to h = 0. 

In the absence of spontaneous symmetry breaking, i.e. for finite systems, and in the absence of external fields, the 
ground state belongs to an (N + l)-dimensional manifold, and can be expanded in the over-complete set of factorized 
states 

\®) = jp I M I '<ty/(M)l*(M)>, (13) 

where / (0, <p) is a weight function. The inner product between states pointing in different directions reads 

(<f>(6',(p') \<$>(0,(t>)) = [cos^cos^ + e^^Osin^sm^l . (14) 



Then, only in the thermodynamic limit we have a set of orthogonal states. Given the continuous SU(2) symmetry 
of the model, spontaneous symmetry breaking implies that the system will select one direction out of all the possible 
choices in the (6*, <fi) space. 

As we are interested in studying the problem of vanishing of peculiar quantum properties, we shall choose initial 
states with given symmetry properties, which, in the finite size limit, do exhibit those properties. 



5 



IV. TRANSITION TO THE SCHRODINGER CAT REGIME 



According to the superposition principle, every linear combination of quantum states is allowed. On the other hand, 
it is well known that superposition cannot be observed on the macroscopic scale because of superselection, the most 
convincing argument being the Schrodinger cat paradox. Then, on this scale, all but a small set of states belonging 
the total Hilbert space are actually forbidden. This process, which leads to a diagonal form of the density operator 
in a preferred basis, implies the vanishing of the most peculiar of quantum properties: state interference. 

In order to study the vanishing of state interference we analyze two different quantities: two-spin entanglement and 
the fluctuation properties of the order parameter M x = (J2 l af) /N. Given a density matrix p, fluctuations statistics 
is associated to the generating function 

G p (A) = Tr{pe^£<^}, (15) 

which is the Fourier transform of the probability distribution function of M x . When the system can be observed in 
m states ^2, • • • ^m, whit related generating functions (A), quantum superposition effects appear if 

AG = Gp (A) VG* (A)^0. (16) 

n— 1 

We expect that in the symmetry broken regime limjy-^oo AG = 0. 

Similar considerations, carried out about entanglement properties, lead to establish that, when N — > oo, only 
factorized states can be observed. 

Even if superselection can be assumed as a principle, the size dependence of quantum interference effects will be 
related to the particular system observed. In the following, we find different decaying behaviors for the XY and the 
XXX models, which are caused by the difference in the symmetry of the two systems. We will start in both models 
by considering symmetric states (which are expected not to survive in the thermodynamic limit) and we shall study 
two-spin entanglement as a function of N. In fact, the existence of degenerate ground state manifolds for any N 
allows to calculate coherence properties as a function of the size of the system. 

For q ubit systems, like spins, two-body entanglement can be measured through concurrence [22j |. As shown in 
Ref. 23], for states which are invariant under the action of the parity operator, the concurrence Cij of two spins at 
sites i and j is related to the quantum correlation functions by simple relations which will be used here. As pointed 
out in Ref. [l9j |. for the XY model at the factorizing field the two-spin concurrence does not depend on the spin 
distance \i — j\. Similar arguments can be used also for the XXX chain. Since we are dealing with superpositions of 
ferromagnetic states, the entanglement will be of ferromagnetic kind as well. In this case, 

Ci 3 ■ = - \pi -P2I -Pin, (17) 
where \pi — p 2 \ is the average value of |ft) (III + III) (ttl) an d Pin is the average value of |tl) (til- 



A. XY model 



Let us first consider the order parameter fluctuations for the symmetric states \a ). These states could be obtained 
by starting with h ^ hp. In this case the exact ground state would have definite parity. For instance, for h > hp, the 
ground state is even. By lowering the field until the value hp is reached, the system is driven in |a + ).The generating 
function is 



G(A,*+) 



G(A,*-)+G (A, 



G(A,* F ,*] 



(18) 



where 



(19) 



It is easy to show that 



G (A, = (^cos A ± i sin A sin 2^ 



iV 



(20) 
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/A \ N 

G(A,*|,^) = (cos — cos20) . (21) 

Then, interference effects (manifested by the off diagonal elements) disappear exponentially with N. 

As a second characterization, we study the concurrence for the symmetric states |a ). This can be easily derived 
using the expression of |a ± ) in terms of j'l'p)- The result (see also Ref. pjl) is given by 

r ( ±\ i o sin 2 2a 

Cij (a ) = (cos2a) —— — w , (22) 

1 ± (cos 2a) 

where the factor (cos29) N ~ 2 derives from the non-orthogonality of j^J) and \^p) and determines the speed of 

classicalization. In the macroscopic limit, C« (or* 1 ) vanishes as [(1 — 7) / (1 +"f)] N ^ 2 - Then, for every finite value of 
the anisotropy 7, entanglement decays exponentially with N. A small but finite anisotropy will enhance entanglement. 
Actually, the 7 = limit implies a non-analytic change 



B. XXX model 



In analogy with the previous case, we introduce a state which is invariant under a given spin rotation. In particular, 
if we choose the state |<I? e ) invariant under rotations about the y axis exp [— iOj^i °f] l^e) = |^e)> we have 



|$e) = i/ d6\* g ), (23) 



r-2n 

where \&g) = ®; |$z (9, 0)), and where \$>i(9,0)) = cos6 |tz) + sin# |4-z). Requiring the normalization of |$ e ) implies 
J dfJ'd^ [cos (0 — 0')] = Af 2 . It is easy to verify that |$ e ) is also an eigenstate of the parity operator. In fact, the 
integration over cancels, in the superposition, all terms with an odd number of down spins. Each would be the 
actual ground state in the presence of an external field directed along the direction 9. 

Let us analyze the order parameter. First, we calculate fluctuations for a given element of the ground state 
degeneracy manifold, obtaining 

1 / A . A 



G® e (A) = j-^ ycos — + i sin — sin 29 j . (24) 

For large N we see that G<s> g (A) ~ exp (iAsin2#). This is the typical expression of the generating function of a non 
fluctuating quantity. Its Fourier transform, which is the probability distribution function of the order parameter is 
indeed, for any 9, a Dirac's delta distribution around sin 29. 
Furthermore, in the superposition state we have 

W = Jp / d6 ' d6 [5 (\ 9, 9')] N , (25) 



where 



g (A, 9, 9') = cos cos (9 -9') + i sin — sin (9 + 9') . (26) 



As N gets large, the vanishing of interference is observed. In the large N limit, the steepest descent method gives 



, A 2 cos 29 , 



(27) 



AG #e (A) = J d9Gs, e (A) 
implying 

AG$ e (A) ~ 1/N (28) 

in the asymptotic regime. 

As far as the two-spin concurrence is concerned, it is straightforward to find 

\ Pl - pa| = / d0de ' [ cos2 6 sin2 6 ' + cos2 °' sin2 9 ] I cos ( e ~ 9')f^ 2 , (29) 
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FIG. 2: Two-spin concurrence as a function of the system size N for the XXX model (dots) and for the XY model (continuous 
lines). In this latter case, we considered two different anisotropy parameters: the red line represents 7 = 0.5, and the blue line 
corresponds to 7 = 0.8. 



and 



yielding 



PIII = / ddd6' cos 9 cos 6 1 sin0 sin & [cos (6 - 6')] N 2 , (30) 



1 / dtidff [tan (0 - 9')} 2 [cos (9 - 9'y N 

2 / d9d9> [cos (9 - 9')} 



C i3 (*.) = " J ^ , (31) 



This result is very simple for N even. In that case one gets 



1 



c % even ^ = 2W^Ty (32) 

This result can be understood taking into account that, given Eq. (TB)) . the inner product between |$(0,O)) and 
|$ (9' , 0)) vanishes exponentially with N . This allows to evaluate integrals, in the large N regime, by means of the 
steepest descent method. It is clear that, when N gets large, [cos (9 — 0')] N is different from zero only for (9 — 9') 
~ 0. Expanding all terms around this value, the concurrence is well approximated by the ratio between two Gaussian 
integrals 

c ^..i ^^K)^ (33) 

which eventually gives Cy (<I> e ) ~ 1/2N in the asymptotic regime. 



V. DISCUSSION 



We tackled the problem of describing how quantum coherence effects vanish as the system size becomes macroscopic. 
Even if this phenomenon is expected to appear in generic systems, two different symmetry-broken model have been 
considered when exact and analytic treatment are possible. In the first one (the XY model in a transverse field), 
because of the discrete symmetry, the ground state spans a two-dimensional manifold. As for the Heisenberg model, 
the dimension of the manifold grows with N, eventually reaching a dense structure. In Fig. [2] we plot the behavior of 
the concurrence Cij for the two models. In the case of the XY chain, we used also different values of the anisotropy 
. A necessary step to determine local quantities, like two-spin entanglement or magnetization, is the introduction of 
the reduced density matrix. Given the peculiar structure of the factorized states we have considered here, calculating 
the reduced density matrix requires the computation of inner products between states defined on (N — m) (where 
?7i is a finite number) spin subspaces) aligned along different directions . Both for the XY and for the isotropic 
model these quantities vanish exponentially in the large N limit, destroying in such a way the quantum interference 
between different components. However, since the XXX model has a continuous symmetry, the ground state manifold 
is continuous as well, and all matrix elements are integrated. The integration implies a reduction in the decoherence 
rate, which turns out to be linear in 1/N. 
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In this paper we considered two exactly solvable models, and studied how superselection tends to destroy their 
quantum properties. A typical tool used to study problems whose solution is not known is the mean-field approxi- 
mation, that, in fact, consists in the introduction of "product states" with the same aspect of those described in this 
paper. For example, in the BCS theory of superconductivity the solution introduced is factorized in the space of the 
modes k. Since in the finite-size case the symmetry is expected to be conserved, while the mean- field states are widely 
unsymmetrical, the linear superposition of degenerate states is a way to restore it. Once the thermodynamic limit is 
performed, all the considerations made in this paper apply. Then, we can conclude that our results apply not only 
to the models explicitly studied, but they could used, in the limit of validity of the mean-field theory, in all systems 
belonging to the classes of symmetry discussed. 
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